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Abstract 

An analysis is developed linking the form of the sound field from a circular source to the radial 
structure of the source, without recourse to far-field or other approximations. It is found that the 
information radiated into the field is limited, with the limit fixed by the wavenumber of source multiplied 
by the source radius (Helmholtz number). The acoustic field is found in terms of the elementary fields 
generated by a set of line sources whose form is given by Chebyshev polynomials of the second kind, and 
whose amplitude is found to be given by weighted integrals of the radial source term. The analysis is 
developed for tonal sources, such as rotors, and, for Helmholtz number less than two, for random disk 
sources. In this case, the analysis yields the cross-spectrum between two points in the acoustic field. The 
analysis is applied to the problems of tonal radiation, random source radiation as a model problem for 
jet noise, and to noise cancellation, as in active control of noise from rotors. It is found that the approach 
gives an accurate model for the radiation problem and explicitly identifies those parts of a source which 
radiate. 

1 Introduction 

A problem in many applications is that of sound generated by circular sources. These include rotors of 
various types such as aircraft propellers and fans, wind turbines and cooling fans; vibrating systems such as 
loudspeakers; ducts such as aircraft engines, ventilation systems and exhausts; and distributed sources with 
axial symmetry such as jets. There have been numerous studies of the noise generation and radiation process 
in each of these areas extending over many decades. These studies can be divided into those which examine 
the relationship between the acoustic source and the physical processes which give rise to it, for example 
the work of Lighthill [1] and of Ffowcs Williams and Hawkings [2] which relate aerodynamic quantities to 
acoustic sources, and those which examine the radiated field generated by a given source distribution, such 
as methods for prediction of the field radiated by pistons and loudspeakers [3-6] or from a known rotating 
source distribution [7-10]. 

There are a number of areas where these issues, those of generation and radiation, overlap. One is the 
general area of source identification. There have been many attempts to develop methods which use acoustic 
measurements to infer, in greater or lesser detail, the source distribution responsible for the acoustic field. 
In the case of rotating sources, some examples include cooling fans [11-13] and propellers [14-17], while a 
number of groups have developed methods for the inverse problem for radiation from a duct termination [18- 
22]. Such studies can have a number of motivations. The first is to use near-field data, for example those 
taken in wind-tunnel tests, to predict the far acoustic field. In this case, the requirement is to extract 
information about source strength and directivity, but there is no need to know which processes generate 
the source. A second motivation, however, is the identification of the noisiest parts of the source with a view 
to reduction of noise at source, for example the identification of "hot spots" caused by unsteady loading on 
a cooling fan [11, 12]. In this case, the link between the aerodynamics and the source is an essential part of 
the solution of the problem. 



1 



In each of the applications of source identification listed, the authors have recognized that the problem is 
(very) ill-conditioned. This can be attributed to physical causes, and is not merely an artifact of the methods 
used. Recent analysis [23-25] has given a framework for the study of this ill-conditioning by quantifying 
the source information which is radiated into the acoustic near and far fields. As described below, it has 
been found that the source can be decomposed into orthogonal modes based on Chebyshev polynomials, 
only a limited number of which radiate a detectable acoustic field, with the limit being fixed by the source 
frequency. 

A second area where the issues of generation and radiation overlap is that of jet noise. Lighthill's acoustic 
analogy [1] is accepted as an exact theory for noise generation by turbulence and there is solid evidence for 
the validity of his source term, as demonstrated by high quality numerical simulation [26]. This knowledge, 
however, is not sufficient to explain certain features of jet noise, in particular the low radiation efficiency of 
subsonic jets and the low order structure of the acoustic field. It is known that subsonic jets radiate only a 
small fraction of the source energy, a view given support by the very small changes in the flow which suffice 
to give large reductions in noise, when control is applied [27]. It is also known that the acoustic far field of 
a jet is significantly simpler than the flow field. In a recent study [28], modal decomposition of the far-field 
noise and of the flow field of a Mach 0.9 jet showed that 24 modes were sufficient to capture 90% of the 
energy of the acoustic field, but 350 were required to resolve 50% of the flow energy. Clearly, a very large 
part of the flow, however energetic it might be, simply does not radiate but it is not obvious if this is due to 
the nature of the source or purely a result of radiation effects. 

The radiation effect has been explained in terms of source cancellation [29, 30] and by viewing the radiation 
process as equivalent to the imposition of a spatial filter using a wavenumber criterion. Such an approach has 
been used by Freund [26] who found that the part of the source which radiates is indeed that part left over 
after applying an appropriate spatial filter. Similarly, Sinayoko and Agarwal [31] apply a linear convolution 
filter to decompose the flow into radiating and non-radiating parts. 

The analysis to be presented below attempts to explain some of these features. Previous work [23-25] has 
found limits on the information radiated from a tonal circular source, motivated by a desire to understand 
the ill-conditioning of source identification methods. These limits have been found without recourse to a 
far-field approximation, making the approach suitable for analysis of general problems. The remainder of 
this paper contains an extension of the theory to explicitly include the radial source term, and to yield 
spectral quantities in the acoustic field of random sources. 

The first extension, which can be viewed as a generalization of previous work on axisymmetric radia- 
tors [32], will help explain radial cancellation effects, which have been studied in jet noise using a far-field 
formulation [30] but not, to the author's knowledge, in the near field. It will be found that for a given 
azimuthal order, many different sources radiate identical acoustic fields, differing only by a scaling factor. 
This result is part of the explanation for the ill-conditioning of identification methods and also opens a 
possible approach to the development of control systems by identifying a class of sources which can give rise 
to practically identical acoustic fields. 

The second extension, to predicting the cross-spectrum between the acoustic pressures radiated by a 
random source to arbitrary points in the near and/or far field, is an extension of an earlier ring-source model 
for radiation from random sources characteristic of jets [33]. In this case, it will be found that the cross- 
spectrum depends on four constants, functions of observer radial separation, which are weighted integrals of 
the source cross-spectrum. 

The results to be presented arise from two different exact theories for radiation from circular sources [24, 
25, 34] which are combined to give a formulation for the information in the acoustic field in terms of radiation 
functions and weighted integrals of the source term. The implications of the results are discussed in terms 
of the information content of the acoustic field and with regard to some of the measurement methods used 
to study noise sources. 
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Figure 1: Coordinate system for disk radiation calculations 
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Figure 2: Transformation to equivalent line source 



2 Tonal disk source 

The problem is initially formulated as that of calculating the acoustic field radiated by a monopole source 
distributed over a circular disk. The system for the analysis is shown in Figure 1 with cylindrical coordinates 
(r, 9, z) for the observer and (a, ip, 0) for the source. All lengths are non-dimensionalized on disk radius. 
The field from one azimuthal mode of the acoustic source, specified as s„(a) expj[m/> — cut], has the form 
P n (k, r, z) expj[n# — cot], with P n given by the Rayleigh integral [23, 35]: 

P n (k,r,z) = J J — ^—^ — dips n (a)ada, (1) 
R' = [r 2 + a 2 - Ira cos V + z 2 ] 1/2 , 
where k is non-dimensional wavenumber (Helmholtz number). 

2.1 Equivalent line source expansion 

The analysis of the nature of the sound field from an arbitrary disk source is based on a transformation of 
the disk to an exactly equivalent line source, an approach which has been used to study transient radiation 
from pistons [3,4], rotor noise [9, 10] and source identification methods [23-25]. 

The transformation to a line source is shown in Figure 2, which shows the new coordinate system (r 2 , 62, z) 
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centred on a sideline of constant radius r. Under this transformation: 



f-r+l Jfcfl' 

P n (k,r,z)= I — K(r,r 2 )r 2 dr 2 , (2) 



2 , „2\V2 



# = {r 2 2 +z 



K ^ = ^j ef) ^s n {a)d6 2 , (3) 
for observer positions with r > 1, with the limits of integration given by: 

^) = cos - il ^ r i. ( 4) 

2 2rr 2 V ; 

Functions of the form of K(r, r 2 ) have been analyzed in previous work [10] and can be written: 

oo 

K(r,r 2 ) = J2^(r)U q (s)(l-s 2 ) 1 / 2 , (5) 

9=0 

where U q (s) is a Chebyshev polynomial of the second kind, s = r 2 —r and the coefficients u q (r) are functions 
of r but not of z. Inserting Equation 5 into Equation 2: 

oo 

P n (k,r,z) = ^2u q (r)C q (k,r,z), (6) 

£ q (k,r,z) = J ^— Uq (s)(r + s)(l-s 2 y/ 2 ds, (7) 

R'= [(r + sf + z 2 f 2 . (8) 

The radiation properties of the integral of Equation 7 have been examined in some detail elsewhere [24, 25] , 
giving an exact result for the in-plane case z — 0: 

£ 9 (fc,r,0)=f(q + l)7rc^^^. (9) 

For large order q, the Bessel function J q (k) is exponentially small for k < q so that the line source modes with 
order q > k generate noise fields of exponentially small amplitude. Since the integrals have their maximum 
in the plane z = 0, Equation 9 says that the whole field is of exponentially small amplitude. This gives an 
indication of how much of a given source distribution radiates into the acoustic field, near or far. 

In previous analyses, two approximations to C q have been developed. One is an asymptotic formula valid 
in the limit k — > oo, derived using the method of stationary phase [24,25]. This will not be required here, 
but we will make use of the far-field form of Equation 7: 



e> kH q + 1 



q R k sin ( 



r + ''/fsTn7 ^ J q +\{ksu\<j)) - jJ 9 (fcsin ( 



(10) 



where R = [r 2 + z 2 ] 1 / 2 and (f> — cos -1 z/R. 

Given the basic information about the form of the radiated field, there remains to establish the relationship 
between the radial structure of the source s„(a) and the line source coefficients u q {r). 
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2.2 Series expansion for spinning sound fields 

A recently derived series [34] for the field radiated by a ring source of radius a can be used to find a second 
expression for the sound radiated by a disk source with arbitrary radial variation: 

i-2tt ^(kR'+nifi) 

Rn= L -^E^ d ^ 



■ 2n+i* 1 y f , m (2n + 4m+l)(2m-l)!! 
4(ai?)V2£-/ > (2n + 2m)!! 

1/2 



X ^+2m+l/2( fc - R )- P n+2m(cOS</))J„ + 2 m+ l/2(fca) J C 11 ) 



with H^\x) the Hankel function of the first kind of order v, J v the Bessel function of the first kind and P™ 
the associated Legendre function. The observer position is specified in the spherical polar coordinates used 
in Equation 10. 

Multiplication by the radial source term as n (a) and integration gives an expression for the field radiated 
by a general source of unit radius and azimuthal order n: 

P (k r z)-\ 2n+1 - f f nm (2rc + 4m+l)(2m-l)!! 

m=0 ^ 

S n+2m (k,r,z) = s n (a)J n+2m+1/2 (ka)H i n 1 l 2m+1/2 (kR)^ ' da. 
Setting z = (<j> = tt/2, P = r): 



Pn(fc, r, 0) = - > -J2 A m S n+ 2m, (12) 
m=0 

1 (2n + 4m + l)(2n + 2m- l)!!(2m- 1)!! 



to! 2 m (2n + 2TO)!! 

where use has been made of the expression [36] : 

pn (M g (2n + 2m-l)!! n - 

2.3 Line source coefficients 

The expressions for P„ from section 2.1 and section 2.2 are both exact and can be equated to derive a system 
of equations relating the coefficients u q (r) to the weighted integrals of the radial source distribution s n {a): 

OO OO j ,, n 

{ £ A mSn + 2m = £ U q (r)j« (« + 1) . (14) 

m=0 g=0 

Under repeated differentiation, Equation 14 becomes a lower triangular system of linear equations which 
connects the coefficients u q (r) and S n+ 2 m '- 



^A m S { :l 2m = J2u g (r)f( q + l) 

m=0 g=0 



e jfer^<?+l(^) 



(v) 

(15) 



k 

where superscript (v) denotes the vth partial derivative with respect to k, evaluated at k — 0. 
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Using standard series [36], the products of special functions can be written: 

e^^^ = ^f>)'+^ t >), (16) 



k j« 



[t/2] 



E U r ) = ^tt E 



21+ 1 j-^ 4 s sl(s + q+ l)\(t - 2s)! ' 



where [t/2] is the largest integer less than or equal to t/2, and 

/a\V2 m /rx 2™+i ~ k 2t+2n+1 ( r 2 \ t 

(-) H^ 1/2 (kr)J n+1/2 (ka) = ( 2 ) E^T- -T K < t(a/r) 



oo , 2t / 2 \ 1 

-(-^Ett -t w ™,*( a A), ( 17 ) 



t! V 4 



with the polynomials V^t and W n ,t given by: 

* /A a ,2s+n+l 

- 5 g r(„ + a + 3/2)r(t- S + n + 3/2)' (18a) 

(18b) 



w»,«(*)=e(! 



„2s+n+l 



_ () v s/ r(n + s + 3/2)r(i - s - n + 1/2) ' 
Given the power series, the derivatives at fc = are readily found: 





fe=o l 







0, w < q; 

j v v\E v - qtq (r), v>q. 



(19a) 



dk? 



{a/rf^H^ (kr)J n+1/2 (ka) 



fe=0 



'0, v = 2v' + 1, u' < n; 

UJ (,pn)!^'-" (a/r) ' « = 2«' + l, ✓ (1Qb) 

Setting u = 0, 1, . . . yields an infinite lower triangular system of equations for u q (r): 

EU = B, (20) 
with U = [uq mi . . .] T and the elements of matrix E and vector B given by: 

B v = I [ T v (r,a)s n {a)da (21b) 
4 Jo 

where 

u = 2i/ + l, v'<n + 2m; 

E Am { fo'-n-2m)! " = 2 ^ + « > - + 2™; (22) 



m=0 



7yWn+2m, t ,'( a / r ) u = 2i/. 



V'! 
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Given a radial source term s n (a), Equation 20 can be solved to find the coefficients u q (r) of the equivalent 
line source modes. Since it is lower triangular, the first few values of u q can be reliably estimated, although 
ill-conditioning prevents accurate solution for arbitrary large q. 

2.4 Radiated field 

From the relationship between the radial source term and the line source coefficients, some general properties 
of the acoustic field can be stated. The first result, already shown in previous work [24, 25] is that, since the 
line source modes with q + 1 > k generate exponentially small fields, the acoustic field has no more than k 
degrees of freedom, in the sense that the radiated field is given by a weighted sum of the fields due to no 
more than k elementary sources. From Equation 20, this result can be extended. 

The first extension comes from the fact that B 2v +i = 0, for v' < n, on the right hand side of Equation 20. 
This means that u q , q = 2v' + 1, is uniquely defined by the lower order coefficients with q < 2v' . The result 
is that the acoustic field of azimuthal order n, whatever might be its radial structure, has no more than k — n 
degrees of freedom, whether in the near or far field. 

A second extension comes from examination of Equation 20. The first few entries of the system of 
equations are: 

"1/2 
r/2 1/4 

: : : 

resulting in the solution: 

w = 2B ; Ul = -2ru = -4rS , (24) 

so that the ratio of Uq and U\ is constant, for arbitrary s n (a). This means that low frequency sources of the 
same radius and azimuthal order generate fields which vary only by a scaling factor, since the higher order 
terms are exponentially small. Again, this result holds in the near and in the far field. 

Finally, if we attempt to isolate a source s n (a) associated with a single line source mode, by setting 
u q = 1 for some q, with all other u q = 0, we find that the line modes must occur in pairs, since if U2 V ' = 1, 
u 2v'+i 7^ 0, being fixed by the condition B 2v >+i = 0, further reducing the number of degrees of freedom or, 
alternatively, worsening the conditioning of the inverse problem. 

2.5 Comparison to far-field methods 

An alternative analysis which is widely used in radiation prediction uses the far field approximations R' w 
R — a sin <p cos ip, 1/R' w 1/R. On this approximation: 

p n ~ (~i) n ~2~^~ J Jn(kasin(t))s n (a)ada, (25) 

so that the radiated field is given by a Hankel transform of the radial source, with a dependence on the polar 
angle <\>. In some sense, this can also be viewed as fixing a limit on the radiated information as in, for example, 
the use of ring sources to study coherence effects on jet noise [30, 33], or as a spatial filter. The approach 
suffers, however, from its inability to give information on the structure of the near field which might be of 
use in understanding such experimental methods as near-field to far-field correlations [37]. The approach 
presented in this paper gives the radiated field, near and far, as the sum of products of two integrals. The 
first of these integrals C q contains only radiation effects while the second u q depends only on the source. 
The source and radiation terms are thus 'uncoupled', simplifying the problem of analysing the radiated field, 
without needing to make a far-field approximation. 
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3 Random disk source 

The second problem considered is that of the noise radiated by a random disk source. This is a general 
problem for broadband noise from rotating systems and is also a model problem for jet noise, extending 
the random ring source problem which has been studied previously in order to examine the effects of source 
coherence on jet noise [33]. The assumptions made are that the source terms are statistically stationary and 
that the statistical properties of the source are symmetric about the source axis. It will also be assumed 
that the non-dimensional wavenumber k < 2, which is a reasonable assumption for the frequency range of 
maximum noise level for a subsonic jet. The result derived is an expression for the cross-spectrum between 
the pressure at two points, which reduces to the power spectrum when the points coincide. The expression 
is quite general and, unlike previous formulae, does not require that the points be in the acoustic far field of 
the source. 

The starting point is an expression for the pressure radiated from a source distributed over a unit disk: 



/•' f' 2 

p(r,6,z,t) = HK ' ' atyda, (2(») 



r2n q(a,ip,t- R/c) 



from which the correlation between p measured at two points (ri, #i, Z\) and (r 2 ,6 2 , z 2 ) is: 



p(n,6 1 ,zi,t)p(r 2 ,6 2 ,z 2 ,t + t) = 

W?L Jo I Jo oio^da^do,. (27) 

Fourier transforming to find the cross spectrum between the points: 

,1 r\ J.2TT ^k(R 2 ~R 1 ) 



W 12 (f) = y—2 / / / — e~5 Qi2(ai,^i;a2,^2)oia2d^idaid^2do2, 

( 47r J Jo Jo Jo Jo R\Ri 



(28) 



/oo 
q(a 1 ,tp 1 ,t)q(a 2 , tp 2 ,t + r)c j27r/r dr, (29) 
-oo 

where Q 12 is the correlation between the source at two points {a\,ipi) and (02,^2), assumed real. 

On the assumption of axial symmetry, the source correlation can depend only on the angular separation 
between two points ip 2 — ipi, so that Q\ 2 and W\ 2 can be expanded in Fourier series in azimuth: 

00 

Qi 2 (oi,Vi;o2,^2)= Yl Q { u\aua2)j m ^-^\ 

m— — oo 

00 

W 12 (r 1 ,e 1 ,z 1 ;r2,e 2 ,Z2)= ]T W&\n, Zl ; r 2 , z 2 y m ^-^\ 

m— — 00 

with: 



Wfr\r 1 ,e 1 ,z 1 ;r 2 ,62,z 2 ) = j^ 2 - f f 

( 47r ) Jo Jo 



•1 r 2v e -j(fc«i+mV'i) 



Ri 



(30) 



1 r 2TT 



e )(kR 2 +mif> 2 ) 



10 JO ^2 

Transforming to the equivalent line source form, as above: 

1 f 1 f 27T C )( kR 2+rn-4> 2 ) 

Q \ 2 ' [ai, a 2 )a 2 d7/>2 da 2 = 

92=0 



Q\ 2 '(ai, a 2 )a 2 dip 2 da 2 



a\ dipt da\ 



— / - Q)% > (a 1 ,a2)a 2 dip 2 da 2 = > u q2 (r 2 ,a 1 )£ q2 (k 7 r2, z 2 ) 

47r Jo Jo K 2 „ _ n 
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which results in: 



W, 



(m) 
12 



1 r 27r p-j(fcfli+m^i) 



Y U H ( ri) r2 ) £ 92 r 2, 22)^ (fc, n, Zi) 



w g2 (r 2 ,ai)ai daid^i, 



where * denotes complex conjugation. The coefficients u qi are found by treating u q2 as the radial source in 
the (ai, - 0i) integral. Up to this point, the analysis is exact but to simplify the development, we introduce 
the assumption k < 2 so that only modes of order and 1 contribute to the acoustic field. 
Solving Equation 20 yields: 

m = 2Bo, 

ui = 4(Sx - rB ), 



with: 



0=1 s m (a)w m (a/r)da, 
Jo 

(a)v m (a) da, 

Jo 



where: 



w m (x) 

Vm(x) 



l_ ^ 1 (2ra + 2 g -l)!!(2 g -l)!! 2q+1 
2ir^q\ 2«(2m + 2 g )!! 



q=0 

x/2n, to = 0; 
0, m^O. 

The result is that the rath azimuthal component of the cross-spectrum between two field points for k < 2 
is given by: 



W[™ ] = C (k,r 2 ,z 2 ) [u < } Cl(k,r 1 ,zi) + u 01 C{ (k, n, Zi)] 
+ d(k,r 2 ,z 2 ) [u w Cl(k,n,Zi) + u 11 £* 1 (k,r 1 ,z 1 )} , 



where 



woo 



u i = 



1 fi 

4 / / Q^\ai,a2)w m (a 1 /r 1 )w m (a2/r 2 )da 1 da2, 
Jo 
1 pi 



Q12 {auO, 2 )w m {a 2 /r 2 )[v m {a 1 ) - r 1 w m (ai/r 1 )] da x da 2 , 



Jo 

1 r l 



u w =8 / Q^\a 1 ,a 2 )w m {a\/ri)[v m {a 2 ) - r 2 w m {a 2 /r 2 )]daida 2 , 
Jo Jo 



(31) 

(32a) 
(32b) 
(32c) 



"11 = 16/ / Q { i2 ) {ai,a 2 )[v m {ai) - r 1 w m {a 1 /r 1 )][v m (a2) - r 2 w m (a 2 /r 2 )] dai da 2 . (32d) 
Jo Jo 

The modal coefficients of the cross-spectrum of a jet noise field, at the wavenumbers of interest in practice, 
are thus fixed by four coefficients, functions of the radial separations r\ and r 2 , which are weighted integrals 
of the source cross-spectrum. 



4 Results 

To check the analyses presented above, and to consider their implications, some results are presented for 
tonal and random disk sources. 



9 



0.1 




5 10 15 
Q 

n = 2 

0.02^ 




5 10 15 



n= 16 

Figure 3: Line source mode coefficients computed using the method of section 2.3 (solid lines) and directly 
from analytical formulae (symbols) for r = 5/4, s = rj, 7 = (circles), 7 = 2 (squares) and 7 = 4 (diamonds) 
for n — 2 and 16. 

4.1 Line source coefficient evaluation 

The first results are a check on the calculation of the coefficients u q (r) comparing those computed using 
Equation 20 and those computed directly from exact closed-form expressions [10] for K{r, T2) in the case 
when the radial source term is a monomial in radius s n = rj. Figure 3 compares the two sets of coefficients 
for 7 = 0,2,4, with the plots terminated at a value of q where the difference between the two sets of results 
becomes noticeable, q w 20. This gives an indication of the effect of the ill conditioning of Equation 20. 
For q < 20, the computed values of u q are reliable. It is noteworthy that for small q, the coefficients are 
practically equal for all values of 7 so that for low frequency radiation, the radiated fields will be practically 
indistinguishable. 

4.2 Tonal radiation from a disk 

As a test of the ability to predict radiation from tonal sources, we present data for the acoustic field of a disk 
source with n = 8, s n = J n (a n ia), where a n \ is the first non-zero root of J n {x). Full numerical integration 
and line source calculations have been performed for two wavenumbers, k = 5 and k = 9, respectively. The 
first 11 line source modes were used in each case, with the modal coefficients being found from Equation 20. 
Sample results are shown in Figure 4, with the data scaled on the value at z = 0, and it is clear that the line 
source model gives accurate results, even when only a subset of the modes is used. From these, and other, 
data, the reliability of the model for tonal sources is confirmed. 

4.3 Low frequency random source 

In order to generate data to test the random disk source model, we must assume a form for the source 
correlation. Michalke [33] gives a form suitable for a ring source which meets the symmetry requirements 
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Figure 4: Acoustic field predicted by full numerical integration (lines) and line source summation (symbols) 
for n = 8, r = 5/4. Real part shown solid, imaginary part dashed. Left hand plot: k = 5; right hand plot: 
k = 9. 

laid out above. With the addition of radially varying terms, Michalke's expression can be extended: 

Qi2(ai,ipi]a 2 ,tp2) = q(a 1 )q(a 2 ) exp 

with a being an azimuthal length scale (3 controlling the correlation in radius. Equation 33 can be inter- 
preted as the product of the local source strengths q{a\) and q(a 2 ) with a coherence function, given by the 
exponentials, which is symmetric in source position and has unit value when the source points coincide. 
The azimuthal components of Qi 2 can be found from mathematical tables [33, 36] as: 

I m (l/a 2 ) (34) 

where I m is a modified Bessel function. 

Figure 5 shows sample results for the predicted cross spectrum between pressure at a point r x = 5/4, 
z\ = and r 2 = 5, < z 2 < 8, for a disk source of unit strength. The reference results are the cross-spectra 
found by full numerical integration of Equation 30. The first comparison is with Equation 31 where the 
functions L q have been evaluated by numerical integration. In the second comparison, the functions L q 
have been evaluated using the exact in-plane result, Equation 9, for z\ = 0, and the far-field approximation, 
Equation 10, for r 2 — 5, < z 2 < 8. All data have been scaled on the numerically evaluated cross-spectrum 
at z 2 = 0. 

The first obvious point from Figure 5 is the similarity of the cross-spectra, even for quite large variations 
in the parameter /?: changing m changes the form of the radiated field, as might be expected, but changes 
in the source correlation have little effect on the radiated field. The second point is that the line source 
approach gives very good results, even for k — 2 where, in principle, the approximation used should start 
to break down. Finally, although computational efficiency is not the primary aim of the method, we note 
that the line source approach converts the four dimensional integral, Equation 30, required at each field 
point, into four two-dimensional integrals which are functions of radial separation only, Equation 32, and 
four one-dimensional integrals Li, giving a large saving in calculation time. 

4.4 Noise cancellation by an equivalent source 

One implication of the results of this paper is that it is not possible to tell different sources apart if, to 
within a scaling factor, they have same line source coefficients u q , for those line source modes with q < k. 
Even without considering errors from background noise or other causes, this is equivalent to a condition on 
weighted integrals of the radial source s n . Any sources which yield the same, or nearly the same, integrals 
B v for v < V, with V a positive integer, in Equation 21, will have indistinguishable acoustic fields for k < V. 
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Figure 5: Cross-spectrum W^™ (-22) scaled on Wf™ (0), ri = 5/4, zi = 0, r 2 = 5. Numerical evaluation 
shown as solid line (real part) and dashed line (imaginary part); Equation 31 with numerical evaluation of 
L q shown as circles; Equation 31 with far-field approximation shown as squares. Parameters: a: k = 1, 
to = 0, a = 1, /? = 100; b: k = 1, to = 0, a = 3, P = 0.01; c: k = 2, m = 1, a = 1, P = 100; d: k = 2, 
m=l, a = 3, P = 0.01. 
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Figure 6: Cancellation effects for radial source terms with n = 2, k = 1, r = 5/4: top figure radiated field 
from original s n {a) (solid) and modified source s n (a) — (s' n (a) (dashed); bottom row source terms s n (a) 
(solid) and (s' n (a) dashed. 
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This conclusion can also be read as a statement about noise cancellation, such as in active noise control. 
The acoustic field of a given source can be cancelled by any source which has the same set of line source 
coefficients. 

An example of this cancellation is shown in Figure 6. The original field is generated using a source term 
s„(o) and the line source coefficients u q of s n are calculated. A secondary source term s' n (a) is generated and 
its line source coefficients u' q are computed. The secondary source s' n is then scaled by a factor ( = uo/u' . 
As a test, s n = J n (a n 2a), with a„2 the second extremum of J n (x), and s' n = 1. The first plot in Figure 6 
shows the field due to s n and that radiated by s n — (s' n . The large reduction, 20dB, in the radiation near the 
source plane is obvious, although there is a small increase in the noise field around z = 1. The source terms 
are shown in the second plot of Figure 6. The secondary source is of much smaller amplitude than s n 
even though it generates a nearly-equivalent field: the effect of matching the line source coefficients has been 
to produce a field which is very similar to that of the original source, even though the source distributions 
are quite different in form and in amplitude. 

5 Conclusions 

The radiation properties of disk sources of arbitrary radial variation have been analyzed to establish the part 
of the source which radiates into the acoustic field, without recourse to a far field approximation. Limits 
have been established on the number of degrees of freedom of the part of the source which radiates and 
the implications of these limits have been discussed for the problems of rotor noise and studies of source 
mechanisms in jets. The analysis has been developed for tonal and for random sources, with implications 
for applications in active control of noise from rotors and experimental analysis of jet noise sources. 
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